Anisotropic fluctuations and quasiparticle excitations in FeSeo.sTeo.s 
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We present data for the temperature dependence of the magnetic penetration depth A(r), heat 
capacity C{T), resistivity p{T) and magnetic torque r for highly homogeneous single crystal samples 
of Fei.oSeo. 44(4) Teg. 56(4). ^{T) was measured down to 200 mK in zero field. We find A(T) follows 
a power law AA '^ T" with n = 2.2 ± 0.1. This is similar to some 122 iron- arsenides and likely 
results from a sign-changing pairing state combined with strong scattering. Magnetic fields of up 
to B =55 T or 14 T were used for the t{B) and C{T)/p{T) measurements respectively. The specific 
heat, resistivity and torque measurements were used to map out the {H,T) phase diagram in this 
material. All three measurements were conducted on exactly the same single crystal sample so 
that the different information revealed by these probes is clearly distinguished. Heat capacity data 
strongly resemble those found for the high Tc cuprates, where strong fluctuation effects wipe-out 
the phase transition at Hc2- Unusually, here we find the fiuctuation effects appear to be strongly 
anisotropic. 



I. INTRODUCTION 

The superconducting iron-chalcogenide compounds 
Fei+ySei-xTe^; have attracted much interest because of 
their many similarities to the high Tc iron-pnictides. The 
transition temperature can be varied between 8K and 
14K by changing the chalcogenide ratio x and reaches 
a maximum at around x — 0.5 {y ~ 0)»i These mate- 
rials share the structural motif of square planar sheets 
of tetrahedrally coordinated Fe with the iron-pnictide 
1111 and 122 compounds (e.g., SmFeAsOi-^Fa; with 
rc=55K and^ and Bai_:rK::,Fe2As2 with Tc up to 38 K^). 
The iron-chalcogenide compounds are structurally sim- 
pler than their iron-pnictide counterparts because there 
are no guest ions or interleaved layers separating the iron- 
chalcogenide layers. Like the iron-pnictides the Fermi 
surface of these compounds is composed of quasi-two- 
dimensional electron and hole pockets located at the 
Brillouin zone corner (M point) and center (F) point 
respectively.— 1^ It is likely therefore, that if the unusual 
Fermi surface topology is the ultimate origin of the pair- 
ing interaction then we might expect the mechanism to 
be similar in these materials. 

In this paper, we report measurements of two aspects 
of the physics of the highest Tc member (x ~ 0.5) of 
the Fei+ySei_2:Te^ series. Firstly, we show data for the 
magnetic penetration depth A(T) in the superconduct- 
ing state. These measurements give information about 
the pairing interactions which drive the superconductiv- 
ity. Secondly, we show measurements of the heat capac- 
ity, magnetic torque and electrical resistivity of exactly 
the same single crystal sample. These measurements are 
used to derive the temperature-field phase diagram of 
this material. By performing the measurements on the 
same sample, the different points on the phase diagram 



measured by these probes is made clear. These mea- 
surements show several unusual features, most notably, 
that unlike some other iron-based superconductors the 
thermodynamics of the superconducting transition are 
dominated by strong thermal fluctuation effects (like the 
high Tc cuprate superconductors). However, unlike other 
superconductors these fluctuations effects appear to be 
strongly anisotropic. 

A powerful way to understand what drives the su- 
perconducting pairing interaction is to determine the k 
dependent structure of the superconducting energy gap 
A. If the pairing interaction v^k' is repulsive in some 
(or all) directions in fe-space, as expected for most spin- 
fiuctuation driven pairing models, then it is likely that A 
will change sign at some point on the Fermi surface." In 
materials where there is a single sheet of Fermi surface 
this inevitably results in gap nodes, i.e., points on the 
Fermi surface where the energy gap is zero. However, 
in the iron-based materials the existence of disconnected 
electron and hole sheet of Fermi surface means that the 
sign-change can be accommodated between two differ- 
ent Fermi surface sheets so that there are no-nodes on 
either^Z^iS. Such as state which is known as the s± state 
appears to be preferred if the scattering is dominated 
by wavevectors close to q = (tt, tt). However if there is 
strong intra-band scattering (with low q) then a state 
with nodes on the hole and / or electron sheets may be 
favored iii"— 

Experiments to determine the gap anisotropy in the 
iron-based compounds have been conducted on many 
different iron-based superconductors. A variety of dif- 
ferent behaviors have been found which might be ex- 
plained by a combination of varying levels of sample dis- 
order and in some cases intrinsic differences in the pair- 
ing state. The first single crystal measurements of the 



temperature dependent penetration depth were reported 
for the 1111 family compomrds PrFeAsOi_^ (Tc=35K>ii 
and SmFeAsOi-a^Fj. (rc=44K).^ These measurements 
showed that for T < 0.2Tc A(T) had a very weak temper- 
ature dependence compatible with a fully gapped state. 
Fitting the full temperature dependence (up to Tc) of the 
superfluid density required a model with two isotropic (s- 
wave) gaps. Experiments on the Co-doped Bal22 com- 
pound, Ba(Fei_2;C0:^)2As2 (Tc =13-24K)l^ showed a ro- 
bust power-law behavior of A(T) oc T" with exponent n 
in the range 2-2.5, which might have been explained by 
a nodal gap in the presence of strong disorder. ^^ How- 
ever, measurements^ of thermal conductivity k did not 
find a substantial value of n/T in the T ~ Q limit as 
would be expected for a sign changing line node in the 
gap and theoretical calculations^^"-^^ suggested that the 
observed power-law in A(r) could result from the in- 
fluence of strong-pair breaking inter-band scattering on 
the intrinsically fully gapped s± state. Strong disor- 
der may be unavoidable in these systems because of the 
doping of Co into the conducting FeAs plane. On the 
other hand, LaFePO is a stoichiometric system in which 
the quasiparticle mean free paths are long enough for 
quantum oscillations to be observable at relatively low 
fieldi^ In LaFePO, A(r) was found to vary close to lin- 
early with temperature which is strongly indicative of 
line nodes in the clean limit i^ Although LaFePO has a 
relatively low T^ of 6K, similar behavior was also seen 
in BaFe2(Aso.66Po.33)2 with rc=30Ki2^ In both cases a 
large value of k/T\t^q was found^^i^ which further sup- 
ports the conclusion that at least some of the Fermi sur- 
face sheets have line nodes in these compounds. So the 
evidence suggests that the pairing state is not the same 
in all iron-based superconductors. It has been suggested 
that the different behaviors are driven by the orbital char- 
acter of the hole sheets which are in turn determined by 
the height of the pnictogen atom.-i^ Kogan^- has recently 
shown that if strong pair-breaking is responsible for the 
T^ power-law behavior of \{T) in the iron-pnictides then 
the coefficient of this T^ term should be related to the 
height of the mean-field part of the specific heat jump at 
Tc and the slope of the upper critical field dHc2/dT. We 
use our data for FeSep.sTeo.s to make a quantitative test 
of this theory. 



II. EXPERIMENTAL DETAILS 

The single crystals used in this work were grown 
from the elements using iodine as the vapor transport 
medium at 680°C targeting the nominal composition 
Fe1.oSeo.5Teo. 5. The details of the synthesis are reported 
elsewherci^ The crystal structure of the exact sample 
measured here (sample 1) was determined by single crys- 
tal x-ray diffraction. No extra iron was observed between 
the Fe(Se/Te) slabs so the mean formula was identified as 
Fei.oSeo.44(4)Teo.56(4). This was additionally confirmed 
by energy dispersive x-ray (EDX) analysis. 



The magnetic penetration depth as function of temper- 
ature has been measured in three different single crystals 
(typical dimensions 300 x 400 x 4/im'^ with the small- 
est dimension along the c-axis), using a radio frequency 
tunnel diode oscillator. ^^ The sample was mounted with 
vacuum grease on a cold finger sapphire rod which was 
inserted into a copper wound solenoid forming the in- 
ductor of a resonant LC circuit with frequency close to 
12 MHz. The RF field was orientated along the c-axis of 
the crystal, thus only in-plane screening currents were ex- 
cited. The experiment was shielded from external fields 
using a mu-metal can and the probe field was < 10~^ T 
thus minimizing any contributions from mobile vortices. 
Two different experimental set-ups were used. One is 
mounted in a pumped ^He cryostat with base temper- 
ature ~1.4K and the second in a dilution refrigerator 
where the sample may be cooled below 100 mK (here the 
resonant frequency is '^14 MHz), ^^ The calibration factor 
relating the measured frequency shifts with temperature 
to AA(T) was determined from the geometry of the sam- 
ple, and the total perturbation to the resonant frequency 
due to the sample, found by withdrawing the sample from 
the coil at low temperature. The procedure is described 
in more detail in Ref. |30. 

The heat capacity of sample 1 from the penetration 
depth study (mass ~ 3/j,g) was measured using a tem- 
perature modulation technique. We use light as a heat 
source to minimize addenda and make the heating uni- 
form over the sample surface, however this does mean 
that the absolute power (and hence heat capacity) are 
difficult to estimate accurately, so our data for C are 
left in arbitrary units. Also the addenda are difficult 
to determine as the length of thermocouple wire which 
contributes depends on the presence or absence of the 
sample. Despite these limitations the method is able 
to measure very small samples of order 1 microgram 
over a wide temperature range with very high resolution. 
These smaller samples are usually much more homo- 
geneous with sharper superconducting transitions than 
larger samples. The sharp transition was critical for the 
accurate characterization of the fluctuation effects as will 
be explained in more detail below. As our present study 
is focused on studying the evolution of the nature of the 
superconducting phase transition with magnetic field the 
lack of absolute scale for C is not a serious limitation. 

The sample was glued with a small amount of VGE- 
7031 varnish to a flattened thermocouple made from 
12/im diameter Chromel and Constantan wires (type 
E).'^^ The sample and thermocouple couple were then 
attached to a temperature controller stage in a vacuum 
enclosure which sits in the ''He bath of a 14 T supercon- 
ducting solenoid. The sample was heated using a mod- 
ulated light source which was generated by an infra-red 
light emitting diode at room temperature and directed 
onto the sample using a plastic optic fibre. The size of 
the temperature oscillations Tac is inversely related to 
the heat capacity C of the samples, Tac = P/^C x /(w), 
where P is the power of the sinusoidal heat source (light) 



at frequency ujM' The function /(w) depends on the 
time constants of the set-up and the operating frequency 
(typically 6-20 Hz) is chosen to be the frequency where 
LoTac has the minimum frequency dependence in the tem- 
perature range of interest. "^^ This ensures that /(w) re- 
mains relatively constant as the experimental time con- 
stants change as a function of temperature and therefore 
C oc T~^ . The temperature of the LED is stabilized to 
minimize drift. The temperature of the stage was moni- 
tored using a Cernox 1050 thermometer, the magnetore- 
sistance of which was corrected for using the data of Ref. 
I33I (although this is negligible in the temperature range 
of interest for the present study). In zero magnetic field 
the thermopower of the thermocouple was taken from the 
standard table for type E thermocouples.-^ The field de- 
pendence of the thermocouple sensitivity was estimated 
by measuring the field dependence of the temperature 
oscillations when a high purity Ag sample (coated with 
a thin layer of black ink to maximize the light absorp- 
tion) replaces the sample on the same thermocouple. At 
T=15 K the sensitivity changes by ~4% between zero and 
14 T whereas at T=24K the change is -2%. 

Electrical resistivity was measured on a small piece cut 
from sample 1 (dimensions 310x 75 x 4/im'^) using the 
standard four probe method with an ac current excitation 
of 100/iA. Contacts were made using pressed Au, which 
resulted in contact resistances of a few ohms. 

Magnetic torque was measured using a piezo resis- 
tive cantilever techniquc'^^ on another small piece cut 
from sample 1 (dimensions 70 x 100 x 4/im'^). Here 
the resistance of an atomic force microscopy cantilever 
is measured using Wheatstone bridge arrangement with 
an ac excitation current of 10/xA at 72 Hz. The sample 
was attached to the cantilever with vacuum grease and 
mounted in vacuum on a rotating platform in the 14 T 
magnet. Checks were made at different current levels 
to ensure that there is no significant self-heating of can- 
tilever and sample in the temperature range of interest. 
Samples from the same batch were also measured using 
the same technique but with pulsed fields up to 55 T in 
Toulouse. Here the sample/cantilever is surrounded by 
^He exchange gas or liquid. 



III. RESULTS 

A. Magnetic penetration depth 

Data for the temperature dependence of A for three of 
our samples of FeSeo.sTco.s are shown in Fig. [T] All three 
samples show a very sharp superconducting transition at 
Tc=13.8K (midpoint) (AT - 0.5 K), and very similar 
low temperature behavior (only sample 1 was measured 
below 1.7K). The absolute temperature dependence of 
the AA(r) is quite reproducible (to within 10%) indi- 
cating that the calibration factor has been accurately de- 
termined and that there are no complications from rough 
edges which can lead to an overestimation of the absolute 
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FIG. 1. (color online) Temperature dependent in-plane pene- 
tration depth for three samples of FeSeo.sTeo.s. For sample 1 
AA is relative to the minimum value measured (at r=0.52 K), 
whereas for samples 2 and 3 the data have been shifted to 
coincide with sample 1 at r=1.7K which was the lowest tem- 
perature to which these samples were measured. The inset 
shows the RF susceptibility (normalized to —1 at T ~ and 
in the normal state) close to T^. The three samples show 
sharp superconducting transitions (ATc ~ 0.5 K). 
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FIG. 2. (color online) Low temperature behavior of the pene- 
tration depth in sample 1 from Fig. [T] The main figure shows 
the temperature dependence of the sample with the back- 
ground from the sapphire rod subtracted (red circles) and 
also with the fitted paramagnetic term subtracted (black di- 
amonds). The lines are fits to a power law (black line) and 
power law plus Curie term (red line) as described in the text. 
The inset shows the measured oscillator frequency shifts for 
the sample plus sapphire rod and also the sapphire rod alone. 



scaleof AA(T),^ 

The low temperature behavior is shown more clearly 
in Fig. [2j Below T ~ 0.5 K there is a small upturn 
in AA(T). Note that the contribution from the sap- 
phire rod on which the sample is mounted is smaller 
and has a significantly weaker temperature dependence 
(see inset to Fig. [5]). The contribution from the rod has 
been subtracted from the data in the main panel. The 
likely origin of this upturn is a Curie-like paramagnetism 
of the normal state. In general the measured penetra- 
tion depth Am is related to the London depth (Al) by 
Am(T) = AL(r)v/l + Xw(T) where xn{T) is the normal 
state susceptibility..^^ Assuming a simple Curie law form 
for Xn{T) then for xn{T) <C 1 the additional contribu- 
tion to the measured penetration depth is given by 
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where rii is the number of magnetic ions per unit cell, Vccii 
is the unit cell volume (=86. 2A'^) and /Lte is the effective 
magnetic moment of the paramagnetic ion. In Fig. [5] we 
show a fit of the low temperature A(T) data to a power 
law dependence plus the AAjvp contribution 
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This equation fits the data very well, with n — 2.2 ± 0.1, 
C = 1.24A/K. This value of C corresponds to an aver- 
age moment of O.lG/is per unit cell (assuming Al(0) — 
5000 A). 38 A fit using a Curie- Weiss form xn =C/{T + 
Tg) gives a slightly better fit (not shown) however as the 
fitted value of Tg = 0. 16 K is below our minimum temper- 
ature this should be viewed with some caution (for this 
fit C = 3.0 A/K and the exponent n was unchanged). 
This small moment could come from the bulk Fe atoms 
or possibly from a small amount of free Fe impurities (al- 
though we note that none were detected in this sample 
by the XRD or EDX measurements). 

The observed power law is close to 2, and as discussed 
above could either indicate a nodal state in the presence 
of disorder or instead could result from strong interband 
scattering between two intrinsically fully gapped Fermi 
surface sheets with sign changing gap between the sheets 
(s± state). Thermal conductivity measurements^S on the 
related compound FeScj, {Tc = 8.8 K) show that k/T\t^q 
is less than 4% of its normal-state value suggesting the 
absence of sign-changing nodes at this composition. As 
far as we are aware there have been no k{T) data pub- 
lished for FeSeo.5Teo.5, so for this composition the exis- 
tence or not of nodes is still an open question. 

Previously, penetration depth measurements have been 
performed by the muon-spin rotation/ relaxation tech- 
nique as well as a tunnel diode oscillator (TDO) tech- 
nique similar to that used here. The /j,SR results of Ben- 
dele et aZ4S for x = 0.5 (Tc=14K) were interpreted using 
a two gap s-wave model, however the data show no evi- 
dence for saturation of A~^(T) at low T due to the lack 




10 12 14 16 18 20 
r(K) 



FIG. 3. (color online) Zero field total heat capacity data. Here 
we have set the constant of proportionality between C and 
T~^^ to unity so the size of the applied temperature oscillations 
can be read off the axis. The solid blue curve is the estimated 
non-electronic background coming from the sample phonon 
heat capacity and also the addenda. The inset shows the 
electronic part of C with the background subtracted together 
with the alpha-model fit (blue curve). The green curve shows 
the same fit convoluted with a Gaussian Tc distribution, T^ — 
13.60 K andar, =0.32K 



of data below T = 1.5K {T/T^ ~ 0.1). The /iSR data of 
Biswas et al^ for a similar composition did show satu- 
ration at low temperature and were also fitted by a two 
gap s-wave model. However, the fact that /j,SR measure- 
ments are conducted in the mixed state means that these 
results are not definitive»2ii^^ The TDO measurements 
of Kim et al^ were reported on samples with x = 0.37 
{Tc ^ 14 K). These data show a power law A ^ T" down 
to the lowest temperatures (T ~ 0.5 K) with n = 2.1 
which is consistent with our findings for FeSeo.sTeo.s- 



B. Heat Capacity 

The heat capacity of sample 1 (from the penetration 
depth study) is shown in Fig. [31 In order to make the 
evolution of the sample's electronic heat capacity with 
field clearer we have estimated the background (which 
includes the samples phonon heat capacity and the ad- 
denda) in the following way. We fitted the zero field data 
over the full temperature range to a polynomial (5th or- 
der) to model the non-electronic terms plus the electronic 
heat capacity calculated using a single gap alpha model. 
In the alpha model^ strong-coupling effects are modelled 
by multiplying the weak-coupling s-wave BCS tempera- 
ture dependent energy gap A(r) by a constant a, and 
then the heat capacity is calculated in the standard way. 
The model works because the temperature dependent 




FIG. 4. (color online) Field dependence of the electronic heat 
capacity (i.e., raw data with the polynomial background in 
Fig. [3] subtracted) . In the upper panel the field is perpen- 
dicular to the planes and increases in 1 T intervals from to 
14 T. In the lower panel the field is parallel to the planes and 
increases in 2 T intervals from to 14 T. Before subtraction 
of the background small (< 1% of the total) corrections were 
applied to normalize all data to the same value at 16 K. This 
is necessary to correct for slight drift in the heating power. 
Field sweeps at this fixed temperature showed that C{H) was 
constant (after correction for the thermocouple field depen- 
dence) within resolution. 



strong-coupling corrections to A(T) are less important 
than the zero temperature value of A(T) for the behav- 
ior of C. For the purposes of the current study we can 
simply regard this as an entropy conserving construction 
for estimating the field independent background. Note 
that the assumption of s or d wave behavior for A(r) 
makes insignificant differences in this range of temper- 
ature (T/Tc ^ 0.5) besides changing a4^ The excellent 
fit of this model to the data is shown in Fig. [31 gives 
^o/Tc — 3.2, which might indicate strong coupling (for 
elemental Pb Aq/Tc = 2.25). 

In the subtracted data some broadening of the super- 
conducting transition is evident, which comes from both 
sample inhomogeneity (a distribution of Tc values) and 
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FIG. 5. (color online) Comparison of the field dependent heat 
capacity Ce{H) for B\\c with the 14 T data for B\\ah. 



also strong thermal fluctuation effects (see later). In or- 
der to quantify the sample inhomogeneity effect in Fig. 
[3] we show a fit to the data with the alpha model con- 
voluted with a Gaussian Tc distribution: i.e., C{T) = 
J-oo C(^' Ta)P{Tc)dTc with the probability distribution 

fmiction P{Tc) ^ exp[-(Tc - T^)^ /{2(j'^J]/ J2TTal^ and 

C{T,Tc) fixed by the above non-convoluted alpha model 
fit. As show in the inset to Fig. [3] this fits the data 
weU with T° = 13.60 K and ctt, = 0.32 K. As fluctua- 
tion effect have not been included here (for simplicity) 
this gives an upper limit on the sample inhomogeneity as 
arjT^ ~ 2%. Note that this distribution of T^ in the 
sample has no intrinsic connection to the scattering pos- 
sible responsible for the T^ behavior of the penetration 
depth. 

The field dependence of the electronic heat capacity 
(i.e., raw data with the polynomial background in Fig. [3] 
subtracted) for the same sample measured in both field 
directions {B\\c and B\\ab) is shown in Fig. ID There 
is a clear striking difference in the behaviors in the two 
field directions. For B\\ab the behavior is close to that 
found for conventional type-Il superconductors. The su- 
perconducting anomaly shifts down in temperature with 
the peak slightly decreasing in height (25% in 14 T), and 
the transition temperature width remains almost con- 
stant. However, in the B\\c direction the anomaly is very 
strongly broadened by the field with the onset tempera- 
ture remaining roughly field independent. The anomaly 
height is reduced by ~ 40% in 14 T. This qualitative dif- 
ference is clearly illustrated in Fig. [5] where we compare 
directly the B\\c and the 14 T B\\ab data. The reduction 
in anomaly height for 14 T B\\ab data is comparable to 
that found for 4T B\\c however, the broadening is not 
comparable at any field. 

The behavior for B\\c is similar to that found on high 
Tc cuprate superconductors where the thermodynam- 



ics of the superconducting transition are dominated by 
strong (critical) thermal fluctuation effects^^."™ Unlike 
conventional low- Tc superconductors where the jump in 
C at Tc is almost a perfect example of a classic second 
order phase transition, in YBa2Cu306.9 {Tc — 93 K) the 
anomaly resembles the A anomaly found at the superfluid 
transition of ^He which is well described by the 3D-XY 
critical fluctuation model. ^" The broadening (with little 
or no reduction in onset temperature) and reduction of 
height of the anomaly is commonly observed in high-T^ 
cuprates4^"— There is some controversy regarding the 
best theoretical picture to describe this. Models based 
on finite size scaling within the critical 3D-XY frame- 
work and fluctuation models based on the lowest Landau 
level (LLL) approximation have both been used.—"— In 
principle, it should be expected that the finite size scal- 
ing should work best at low fields (compared to Hc2 at 
T = 0) and the LLL model at high fields;^^ The finite size 
effects arise because the magnetic field introduces a new 
length scale ds = (^o/Tri?)^^^, which limits the growth of 
the physical size of the fluctuations as T approaches Tc- 
When the transition is dominated by such strong fluc- 
tuation effects the true phase transition at Hc2 is wiped 
out as the fluctuations in the order parameter are much 
larger than the mean value over a large range of {H, T) 
space close to the Hc2{T) line.^^ The heat capacity still 
displays a broad hump at the 'i?c2(r)' line but the resis- 
tivity for example does not become zero until lower tem- 
perature when the vortex lattice freezes and phase coher- 
ence is achieved. The strong fluctuation effects in high Tc 
cuprates are driven by their short coherence lengths, low 
dimensionality and high critical temperatures i^ Com- 
pared to the cuprates FeSco.sTeo.s is relatively isotropic 
and has a much lower Tc but does have a comparably 
short coherence length. The estimated irreversibility fleld 
at T = of FeSco.sTeo.s is similar to that found for un- 
derdoped YBazCuaOg.s (Ref. [H) or YBa2Cu408 (Ref. 
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FIG. 6. (color online) Resistivity versus temperature for sam- 
ple 1 for fields between and 14 T applied either parallel or 
perpendicular to c. The inset in the upper panel shows a di- 
rect comparison of the resistivity (left hand scale) and heat 
capacity (right hand scale) of the same sample (sample 1) in 
zero magnetic field. 



C. Resistivity 



The fleld broadening of C{T) for B\\c is in contrast 
to the more conventional behavior observed in 122 com- 
pounds such as BaFe2(Aso.66Po.33)2 with Tc=30K,^^ 
however evidence for strong fluctuation effects have been 
seen in the field dependent heat capacity of the 1111 com- 
pound NdFeAsOi_:rF^ (Ref. |5J). An unusual aspect of 
the behavior of FeSeo.sTeg.s is that the fluctuations ef- 
fects appear to be strongly anisotropic. In YBa2Cu306.9 
although the 'Hc2{Ty and irreversibility fleld iJirr are 
anisotropic, the behaviors in the two directions are iden- 
tical apart from a difference in the fleld scale (the heat 
capacity for IT ||c is identical to 8T||a&)42ii^ Here the 
behaviors are qualitatively different in the two field di- 
rections. Fluctuations effects are minimal for B\\ab but 
strong for B\\c. It would appear that in-plane field does 
not disrupt the phase coherence between the planes but 
for B\\c the effective dimensionality is strongly reduced 
and fluctuations effects enhanced. 



The temperature and fleld dependence of the resistiv- 
ity of sample 1 is shown in Fig. [6l Although there is again 
more broadening of the transition for B\\c compared to 
B\\ab it is much less evident than that seen in the heat 
capacity data. In the inset to Fig. IH] we show a direct 
comparison between the zero fleld heat capacity and re- 
sistivity data for the same sample. The mid point of the 
resistive transition is ~0.75K higher than the mid-point 
of the heat capacity. The peak in the heat capacity is a 
further 0.5 K lower and corresponds to the point where 
the resistivity is ^1% of its value just above Tc- This 
difference is understandable as a large reduction in the 
resistivity will be achieved when a sufficient amount of 
the sample to create a percolation path becomes super- 
conducting. Often this amount will be much less than 
50% but depends on the spatial distribution of Tc in the 
sample. 
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FIG. 7. (color online) Torque hysteresis loops. The left two 
panels show data in quasi-dc field for a small piece cut from 
sample 1. The right two panels show data taken in pulsed 
field for a second sample. The field directions are actually 
~ 5° from the symmetry directions indicated on the figure. 



D. Torque 

Torque data in quasi-static fields up to 14 T for a small 
piece cut from the heat capacity sample is shown in Fig. 
[7]at angles of approximately 5° from B\\c and -B||a&. The 
angle offset is necessary because at the symmetry points 
the torque is exactly zero. These small offset make negli- 
gible difference to the derived irreversibility fields as the 
anisotropy is not strong. Data taken in pulsed fields for a 
different sample from the same batch is shown in Fig. [T] 
In both cases the torque is strongly hysteretic indicating 
strong pinning. In the pulsed field B\\c data oscillations 
are evident which likely arise from fiux jumps. In all data, 
the torque loops close at within our resolution at an ir- 
reversibility field (iJhr), and no further kinks at higher 
field (which might mark Hc2) are evident. 



E. Phase diagram 

A compilation of the data derived from the heat capac- 
ity, resistivity and torque measurements is plotted in Fig. 
|S]as a function of reduced temperature Tc{B)/Tc{B — 0). 
It is important to realize that in general these different 
experimental techniques will measure different points in 
the {H, T) phase diagram. For a superconductor with- 
out strong fiuctuation effects there will be a step in 
the heat capacity at Tc which will be broadened by in- 
homogeneity and will shift down in temperature under 
applied field marking the Hc2(T) line. The width of 
the broadened step in field will be related to the width 
in zero field and the slope of Hc2{T), i.e., ATc{H) — 
Hc2ATc{0)/{dHc2/dT). In this case the mid point of the 
rise of the C{T) can be used to determine Hc2{T). How- 
ever, in the present case, as the broadening of C{T) in 
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FIG. 8. (color online) H ~ T phase diagram for FeSeo.sTeo.s 
close to Tc derived from heat capacity, resistivity and torque 
measurements. For the heat capacity the temperature of 
the peak in C{T) and the mid point of the rise above the 
peak the transition are plotted. For the resistivity, Rso% 
and Ri% are the points at which the resistivity has fallen 
to 50% and 1% of extrapolated normal state resistivity above 
Tc respectively. For the torque the point where the hystere- 
sis loop closes (within our resolution), i.e, the irreversibility 
field fioHirr is plotted. The lines are fits to the power-law 
H = Hoil - T/Tc)", to guide the eye. 



field cannot be explained by inhomogeneity alone there 
must also be strong thermal fluctuation effects as de- 
scribed above. In the case of phase transitions where 
fluctuation effects are strong, the peak in C{T) is usually 
used to deflne the phase transition point, and so for the 
present case this likely gives the best estimate of Hc2{T). 
At the point at which the magnetization becomes irre- 
versible (the irreversibility line) the sample will acquire 
a finite critical current. Hence the irreversibility line as 
measured by our torque measurements marks the point 
where the critical current takes a finite value - the size of 
which depends on the resolution of the torque measure- 
ments. It should be expected that this irreversibility line 
is close to the point where the resistivity becomes zero 



as this is also where the critical current becomes equal 
to the resistivity measurement current. In the case that 
the irreversibility line and Hc2 are not in close proximity 
it is difficult to know what point in the resistivity curve 
corresponds to Hc2{T). For example, if the vortex vis- 
cosity is very low the resistivity in the mixed state will be 
close to the normal state value and there will be almost 
no anomaly at Hc2[T)^^^ 

In the present case, the peak in C{T) coincides reason- 
ably well with the point on the resistivity curve where the 
resistance has fallen to 50% of its extrapolated normal 
state value (i?5o%). This is the case for both field direc- 
tions (see Fig. [8]). Hence is it reasonable to conclude that 
both of these two quantities are measuring the tempera- 
ture dependence of Hc2- The 1% resistivity point (i?i%) 
is better correlated with the torque irreversibility field as 
expected from the above discussion. 

For B\\c there is clearly a very marked difference be- 
tween the irreversibility line derived from the torque mea- 
surements and the thermodynamic transition measured 
by specific heat. The average slope of the torque irre- 
versibility field dHl^JdT = -44T/K whereas for the 
peak in C{T) it is -95T/K and the midpoint -150 T/K. 
Hence estimates of Hc2 from the irreversibility field for 
this field direction will serious underestimate the true 
values. Although the data are better fitted by a power 
law H = {1 — T/Tc)'^ as indicated in the figure rather 
than a linear relation, the average slope still provides a 
useful rough comparison. The large separation between 
the irreversibility field and the peak in the specific heat 
which we expect to be close to the mean-field Hc2 line 
is again reminiscent of the behavior found in the high- 
Tc cuprates and is strongly suggestive of strong fluctua- 
tion effects4^"— For B\\ah the lines are much closer to- 
gether. The corresponding values for the average slopes 
are -90 T/K, -160 T/K and 200 T/K respectively. 

The anisotropy of Hc2 derived from the specific heat 
measurements is considerably less than that derived from 
the torque irreversibility field. Taking the peak in C{T) 
as the most reliable measure of Hc2{T) we see that for 
B\\c it is almost linear with temperature whereas for 
B\\ab there is some curvature close to Tc{Q)- Above 
^ 8 T where hI^ (T) is linear the anisotropy is small. 
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Braithwaite et alM- have also reported a difference be- 
tween critical fields determined by resistivity or heat ca- 
pacity measurements for FeSeo.48Teo.52. In their case, 
they found i?5o% and the mid-point of the C(T) curve ap- 
proximately coincided for H\\ab but not for H\\c. Their 
heat capacity data showed almost no anisotropy. C{T) 
data for fioH=9T was almost independent of the field 
direction. It is likely that the differences in the H\\ab 
data compared to ours is due to the much wider transi- 
tion width of their sample. The zero field C{T) data in 
Ref. [53 are approximately three times wider than those 
here. The field broadened transition width of our data 
for fioH—9T ||c is comparable to the zero field width in 
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FIG. 9. (color online) H — T phase diagram for FeSeo.sTeo.s 
over the full range. The torque data mark the irreversibly 
line whereas the heat capacity data C{T) are two different 
estimates of the behavior of Hc2- The lines are guides to the 
eye. 



Ref. l57|. 

The phase diagram over a wide temperature-field range 
is shown in Fig. IHJ As reported previously^2i^ the irre- 
versibility fields for the two field directions cross at low 
temperature with the anisotropy approaching unity. This 
and the general shape of these H[^^{T) curves, has been 
interpreted as evidence that Hc2 is Pauli limitedi ^^i^^ 
However, given the difference between Hi^^ and /Jc2 re- 
vealed by the lower field heat capacity data, it is pos- 
sible that this high field behavior is also strongly infiu- 
enced by thermal ffuctuation effects and vortex pinning. 
High field heat capacity measurements will be necessary 
to determine whether the temperature dependence and 
anisotropy of Hc2 follow that of iJirr at high field and low 
temperature. 



IV. PAIR-BREAKING 

As mentioned in the introduction, Kogan^i has cal- 
culated that in anisotropic superconductors with strong 
pair-breaking there should be a universal relationship be- 
tween the co-efficient of the T^ term in the low tempera- 
ture penetration depth, the height of the mean-field part 
of the specific heat jump ACmf, and the slope of the Hc2 
near Tj.. In SI units the relation is 



2rT^4 



n = 



ACmf A'T, 

I dHa'. 



dT It, 



90_ 

47r 



(3) 



Note here that ACmf is units of J m~^ K~^ and (po is 
the flux quantum. Fixing n = 2 in Eq. [5] we get A — 
10.3 A/K for the T^ slope of the low temperature A(T) 
data. As argued above, the peak in C{T, B) is likely 



to provide the best estimate of Hc2{T) and from this 
we get A«o^ = 6.9 T/K for B\\c and T between 0.85 
and 0.97 Tc. As our specific heat data is not in absolute 
units we estimate ACmf from the data of Tsurkan et 
al^ for the specific heat of FeSeo.sTco.s. This data is 
very similar to our own with regard to sharpness and 
value of Tc and slope of Hc2, and from this we estimate 
ACmf = 0.42 Jmor^K^^ Taking the mid point of the 
heat capacity transition T — 13.6 K as Tc, we get from 
Eq. [21 47rf7/0o — 0.65 ± 0.15 which is indeed close to 
unity as predicted by the Kogan theory (the uncertainty 
is dominated by that in Tc). Hence, it is likely that the 
T^ behavior of the penetration depth does result from 
strong pair breaking. However, this does not rule out 
the possibility of intrinsic nodes if the material could be 
made sufficiently clean. 



depth data display a power law behavior AA oc T" with 
n = 2.2±0.1, which is similar to some 122 iron-arsenides. 
It is likely this results from the influence of strong pair 
breaking with a sign-changing pairing state which may 
or may not have intrinsic nodes. ^^"^^ Heat capacity mea- 
surements have shown clear evidence for the presence of 
strong thermal fluctuations, at least for relatively low 
fields {^jlqH < 14 T). The behavior is reminiscent of the 
high Tc cuprates, and like in these compounds we show 
that there is a sizable vortex liquid regime between the 
irreversibility line and the upper critical field. Unusually, 
the fluctuations influence the behavior with field along c 
much greater than for field along the planes. 
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